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We have analyzed 44.3M single muons collected by MACRO from 1991 through 2000 in 2,145 live 
days of operation. We have searched for the solar diurnal, apparent sidereal, and pseudo-sidereal 
modulation of the underground muon rate by computing hourly deviations of the muon rate from 
6 month averages. We find evidence for statistically significant modulations with the solar diurnal 
and the sidereal periods. The amplitudes of these modulations are < 0.1%, and are at the limit of 
the detector statistics. The pseudo-sidereal modulation is not statistically significant. 

The solar diurnal modulation is due to the daily atmospheric temperature variations at 20 km, the 
altitude of primary cosmic ray interactions with the atmosphere; MACRO is the deepest experiment 
to report this result. The sidereal modulation is in addition to the expected Compton-Getting 
modulation due to solar system motion relative to the Local Standard of Rest; it represents motion 
of the solar system with respect to the galactic cosmic rays toward the Perseus spiral arm. 

PACS numbers: 95.85.Ry, 98.70.Sa 


I. INTRODUCTION 

MACRO is a deep underground experiment 
designed to search for monoples and rare par¬ 
ticles in the cosmic rays. Because of its great 
depth and sensitivity to high energy primaries, 
the apparatus also probes the physics of cosmic 
rays with energies in excess of 10 TeV. Here we 
describe a search for the first harmonics of the 
sidereal and solar diurnal modulations in the un¬ 
derground muon rate using the MACRO detec¬ 
tor. Sidereal modulations result from the mo¬ 
tion of the solar system with respect to a lo¬ 
cally isotropic population of cosmic rays. For in¬ 
stance, a local source of cosmic rays could intro¬ 
duce a sidereal wave into the underground muon 
rate. Or, if the local density of halo cosmic rays 
were sufficiently high, a sidereal wave would be 
imprinted on the underground cosmic ray rate 
as the solar system moved through them. Solar 
diurnal modulations at energies greater than 10 
TeV on the other hand are mostly the result of 
meteorological effects at or near the primary site 


II. THE FIRST HARMONIC SIDEREAL 
AND SOLAR DIURNAL MODULATIONS 

There are two primary processes that can 
lead to periodic modulations in the underground 
muon rate. One process, the Compton-Getting 
effect l]^, is due to an observer’s motion with 
respect to a locally isotropic population of cos¬ 
mic rays. The modulation is introduced as the 
Earth’s rotation periodically turns an observer 
toward and away from the direction of motion. 
Since the Earth rotates on its axis in a sidereal 
day, the Compton-Getting modulation is seen 
when muons are binned in sidereal time. The 
second process that leads to a periodic mod¬ 
ulation in the muon rate deep underground is 
the result of the competition between pion de¬ 
cay and interaction in the upper atmosphere. 
As the atmosphere cools at night (or during the 
winter), the density increases and the pions pro¬ 
duced in cosmic ray collisions with the atmo¬ 
sphere are fractionally more likely to interact 
than decay than during the day (or summer) 


muon rate. The daily modulation introduced by 
solar heating and cooling is seen when muons are 
binned in solar-diurnal time; yearly or seasonal 
modulations are seen when events are binned 
with a period of a (tropical) year |^, ^ . 

Underground experiments measure the muon 
rate R. So a sensitive search for sidereal and 
solar diurnal modulations in MACRO is a search 
for periodic variations in the underground muon 
rate as compared with the average rate, 

AR/Ro= jj{AI/Io)dndA, (1) 

where AR = R — Rq, R is the measured cos¬ 
mic ray rate, Rq is the average cosmic ray rate, 
AI/Iq = {I — Io)/Iq is the deviation from the 
average cosmic ray intensity, and the integration 
is over the acceptance of the detector. 

A. The Compton-Getting Effect 

An observer moving with velocity v relative to 
the rest frame of a cosmic-ray plasma will detect 
a deviation due to this motion from the average 
cosmic-ray intensity, an effect first described by 
Compton and Getting |^. If the cosmic rays 
have a differential power-law energy spectrum, 
, then the first harmonic of this deviation is 
given by 

AI{iP)/Iq = ( 2 -I- 7 ) (u/c) cos/3 cos 5 

X cos27r[(A — a)/24], (2) 

where AI{ijj) = I{tp) — /q; I{tp) is the cos¬ 
mic ray intensity in a direction with space an¬ 
gle ijj between the direction of the detector’s 
maximum sensitivity, (/3, A), and the direction 
of motion, (<5, a), in celestial coordinates on the 
sky; Iq is the average cosmic ray intensity; and 
(2 -I- 7 )(u/c) is the first harmonic amplitude Q. 
In this expression, /3 is the declination and A is 
the right ascension (in hours) of maximum de¬ 
tector sensitivity, and S is the declination and 
a is the right ascension (in hours) toward the 
direction of motion. Since a cosmic-ray detec¬ 
tor is fixed on the rotating earth, it sweeps the 
sky with period equal to a sidereal day. Under 
these circumstances A = A(<), and the Compton- 
Getting signal is modulated with a period of a 
sidereal day. 

Surface detectors are typically most sensitive 
to the declination crossing the local zenith and 
the right ascension crossing the meridian. For 
such a detector, /3 is equal to the detector’s lat- 


The A response of the detector therefore implies 
that a Compton-Getting signal will be modu¬ 
lated in sidereal time, t*, by the phase term 

27r[(A-a)/24] = 27r(r*-a/24), (3) 

where the sidereal phase t* = t*/24. 

As a deep underground detector with 
angular-dependent acceptance and overburden, 
macro’s directional sensitivity (/3, A) to side¬ 
real signals differs from the ideal surface detec¬ 
tor described above in two ways. First, the de¬ 
tector response may be peaked away from the 
meridian. In this case, MACRO would detect 
the signal at a sidereal hour displaced from its 
true right ascension. The upper panel of Fig. 1 
shows a histogram of the hour angle of the arriv¬ 
ing muons; this histogram shows that the mean 
of the hour angle distribution falls at —7.7° from 
the meridian, an offset that translates to a cor¬ 
rection of —31”’' to the sidereal hour of the de¬ 
tected signal maximum. 

Second, the detector response may be peaked 
away from the zenith. For a detector whose re¬ 
sponse peaks away from the zenith, the decli¬ 
nation of the detector’s maximum sensitivity in 
eq.(||) is not simply its latitude and (3 needs to 
be replaced by an ‘effective declination’, /3e//- 
It is evident from eq.(|^ that this replacement 
affects only the amplitude of the signal and not 
the sidereal hour of its maximum. For MACRO 
Peff can in principle be determined by a peak 
in the declination distribution for the muons 
used in the analysis. The lower panel of Fig. 1 
shows this declination distribution. The situa¬ 
tion is clearly complicated by the broad distri¬ 
bution with mulitple peaks at 15°, 42°, and 64°. 
(macro’s latitude = 42.4528°.) In the analysis 
presented here, we use the mean of the distribu¬ 
tion as the effective declination, (3eff = 37°. It 
should be noted, however, that a value of (3eff 
between 15° — 42° only alters the signal ampli¬ 
tude by less than 30%. 

Possible Compton-Getting signals with a side¬ 
real period might result from solar motion with 
respect to galactic sources of cosmic rays |||] or 
from motion with respect to halo cosmic rays 
residing in the galactic disk. The effect of solar 
system motions on the underground muon rate 
gives the modulated cosmic ray intensity 

[A/(r*)//o]cG = {2 + ^){v/c) cos(3eff cos 5 
X cos27r(T* — a/24). (4) 

Thus the Compton-Getting signal expected in 
the undergound muon rate would have the form 

TAn/ \/7-il ...fA / lc\ A \ / r'\ 


One component of the Compton-Getting signal 
seen by MACRO in sidereal time is due to the 
solar motion relative to the Local Standard of 
Rest (LSR) [|j , the rest frame of the local cosmic 
rays. Adopting v = 16.5 km/s toward I = 53°, 
b = 25° for this motion , we compute = 

2.6 X for the amplitude of this wave and 
a^SR — 17 . 8 ^ for its right ascension. 

Previous searches for this effect have given dif¬ 
ferent results depending on the energy of the 
primaries probed. Underground muon obser¬ 
vatories and air shower arrays (E > 10 TeV) 
find sidereal modulations with amplitudes in the 
range 5 x 10“"^ — 10“^, consistent with a drift 
velocity of a few hundred km/s [Q. Shallow un- 
dergound muon telescopes and neutron monitors 
(E < 500 GeV) observe statistically significant 
sidereal modulations that are explained by solar 
wind effects l^. Underground muons detected 
by MACRO have energies at the surface in ex¬ 
cess of > 1.25 TeV and are due to interac¬ 
tions of primaries whose energy exceeds 10 TeV. 
Therefore, the observed muon rates are unaf¬ 
fected by the solar wind. 

If cosmic ray muons deep underground are 
binned in local solar diurnal time, t©, there will 
be a Compton-Getting signal that is a conse¬ 
quence of the earth’s orbital motion through the 
local solar system cosmic rays. Since the muons 
have energies at the surface > 1.25 TeV that 
leave them unaffected by the solar wind, they 
are locally isotropic when averaged over the so¬ 
lar diurnal period. Ignoring the earth’s orbital 
eccentricity, the Compton-Getting modulation 
due to this orbital motion peaks near 6 am local 
time when the detector is sensitive in a direction 
most nearly parallel to the earth’s orbital veloc¬ 
ity, and the cosmic ray intensity is modulated in 
solar diurnal time by the phase term 

27r(A - ag) = 2 tt{tq - 6/24), (6) 

where the solar diurnal phase Tg = tg/24 Q. 
In this case, 

[^I{tq)/Io]cg = {2 +-f){v/c) cosPeff 

X cos27r(rg — 0.25). (7) 

The Compton-Getting signal in solar diurnal 
time would have the form 

[Ai?(rg)/i?o]cG = ATcg cos 27r(Tg - 0.25) 

= 4.55 X 10“'^ cos/3e// 

X cos27r(rg — 0.25), (8) 

where the first harmonic amplitude has been 
computed from the earth’s average orbital 


B. Atmospheric Effects 

In addition to the Compton-Getting modula¬ 
tion, atmospheric effects modulate the muons in 
solar diurnal time, [AI^^tq)/I o]atm- This atmo¬ 
spheric modulation is the result of density vari¬ 
ations at the altitude of first interaction - when 
cold at night, the density at the primary inter¬ 
action site is fractionally higher than the daily 
average and pions are fractionally more likely to 
interact than decay, thereby giving fewer muons 
relative to the average underground muon rate. 
During the day, the reverse occurs. The mod¬ 
ulation of the underground muon rate in solar 
diurnal time due to atmospheric density varia¬ 
tions is given by 


[Ai?(rg)/i?o ]0 = A:gCOS27r(Tg - ^g), (9) 

where ^g is the phase of maximum modula¬ 
tion. Atmospheric effects also modulate the 
underground muon intensity on a yearly time 
scale. This is the so-called ‘seasonal modula¬ 
tion’, [AI{Tyr)/Io]yr, where Tyr is the yearly 
phase. If Nq is the number of solar days elapsed 
since the beginning of the year and Ngoi is the 
number of days in a tropical year, then Tyr = 
(Vg-|-tg/24)/Vso; and the seasonal modulation 
is proportional to cos 2 t:(T yr — ^yr)- The sea¬ 
sonal modulation has the form 


[AR(Tg )/Ro]yr — Kyr COS27r(Tyr ^yr) • (fH) 

In Fig. 2 we show for MACRO the percent¬ 
age deviations of the monthly average muon 
rate from the yearly average (solid circles), as a 
measure of [Ai?(Tg)/i?^r, for the years 1991- 
1997. The fit of eq.(|l^ to these data gives 
Kyr = 0.011 ± 0.0004 and ^yr = 0.54 ± 0.06. 
Superposed on the muon data are the devia¬ 
tions in the mean monthly ‘effective tempera¬ 
ture’ of the atmosphere from the yearly aver¬ 
age, ATf.jj/ < Teff > (open triangles), as de¬ 
fined in . We have repeated the analysis from 
Q| that computes the correlation of the mean 
monthly deviation in the muon rate with the 
mean monthly deviation in the effective temper¬ 
ature of the atmosphere, 


[Ai?(rg)/i?o]yr = ttrATeff/ < Teff > . (11) 

A regression analysis shows that ay = 0.91 ± 
0.07, or equivalently 0.42±0.03 %/K, consistent 

•<i 1 .J ..r. 1 * 


C. The Pseudo-Sidereal Modulation 


As first pointed out by Compton and Get¬ 
ting Q, atmospheric effects can mimic a side¬ 
real modulation if the yearly (seasonal) wave 
beats with a solar diurnal wave. Under these 
circumstances, the beats due to atmospheric ef¬ 
fects have the form 

[2Ai?/i?o]beats = KQKyr X 

COS27r[(A^so/ \^Tyf ^0 Cyr] 

-\-KQKyT X 

COS27r[(A^5o/ l^Tyr ^0 -f ^yr\ 

— ^Q^yr X 

COs27r(r* - - Cyr) 

-\-KQKyr X 

cos 27r(Tpse?ido Co d“ Cyr)- (f^) 

In this expression, the first modulation, with 
a period of {Ngoi + 1) solar days, represents a 
pseudo-sidereal modulation resulting from at¬ 
mospheric effects. This modulation was named 
the ‘antisidereal’ modulation by Farley and 
Storey and is often called the ‘antiside¬ 
real wave’ in the literature. The second mod¬ 
ulation, with a period of {Ngoi — 1), solar days 
represents a modulation in unphysical ‘pseudo- 
sidereal time, ipseudo ^sol ^ ^q/^^^soI f)j 
and is purely the result of the beating of the 
yearly and daily solar modulations. 

To correct for the pseudo-sidereal effect, the 
method described by Farley and Storey |j^ is 
often used. In this method, the amplitude of 
the pseudo-sidereal wave, which is equal to the 
amplitude of the pseudo-sidereal wave, is deter¬ 
mined by binning the data in pseudo-sidereal 
time. The phase of the pseudo-sidereal wave 
is determined by reflecting the pseudo-sidereal 
wave vector about the solar wave vector. 


III. DATA ANALYSIS 

This analysis presented here requires only the 
arrival time of each muon. In principle, the anal¬ 
ysis is simple: the arrival time of each muon is 
binned separately in histograms of sidereal, solar 
diurnal, and pseudo-sidereal time; an accurate 
account of the live-time in each histogram bin is 
accumulated; at the end, the number of events 
in each bin is divided by the live-time in that 
bin to give a muon rate. The results are based 
upon separate analyses of the muon rate in each 
histogram bin compared with the average rate 


For the sidereal analysis, it must be empha¬ 
sized that the results are not based on tracking 
information, that is, the determination of the 
right ascension and declination of each muon. 
Here is the aim of the sidereal analysis. If 
the muon rate were completely uniform over 
the celestial sphere, the detector would mea¬ 
sure a uniform rate underground as a function 
of sidereal hour. If there were an excess of 
muons from some particular celestial direction, 
MACRO would measure an increase in the muon 
rate as the detector’s maximum response swept 
over that direction. In the sidereal analysis, we 
are searching for this excess. 

A. Run/Event selection 

In this analysis we included data runs start¬ 
ing from the beginning of MACRO data taking 
with 6 supermodules in November 1991 through 
May 2000. (A complete description of the de¬ 
tector during its running is given in . ) A typ¬ 
ical data run had a duration of 6-12 hours, and 
over 9,500 data runs were included in this anal¬ 
ysis. The analysis proceeded by first dividing 
these individual data runs into 17 run sets of 
approximately 6 months duration during which 
the detector acceptance remained constant. We 
then filtered the data to include only muons with 
single tracks in both views. In this analysis, 
we accept events that fire at least 3 planes in 
the lower detector. For events crossing only 3 
planes, the trigger efficiency is approximately 
72%; for events crossing more than 3 planes, 
the trigger efficiency is > 99%. Once filtered, 
we compiled for each run set a histogram of the 
single muon rate per hour for all runs in that run 
set and then fitted a Gaussian to the resulting 
distribution. 

We implemented run cuts as follows. A run 
was excluded from further analysis if: 

• not all 6 supermodules were active, or 

• the wire efficiencies were < 90% and/or 
the strip efficiencies were < 80% during 
the run, or 

• the single muon rate for the run was > 5cr 
from the mean single muon rate for that 
run set. 

For the efficiency cuts, wire and strip efficien¬ 
cies were determined for each run by computing 
the average number of wires and strips recording 
hits for all single muons crossing 10 planes. 

There are 44.3M muons in the total data set 


eluded runs is 2,145 days. The average sin¬ 
gle muon rate in the total data set is 860.53 
muons/hour /6 supermodules. 


B. Histograms of deviations from the 
mean solar diurnal, apparent sidereal, and 
pseudo-sidereal muon rates 

We searched for the solar diurnal, apparent 
sidereal, and pseudo-sidereal modulations as fol¬ 
lows. First, event histograms for each run were 
constructed for the three periods by binning 
the arrival time of each muon according to its: 

(1) local solar diurnal time at the Gran Sasso; 

(2) local sidereal time; and (3) local pseudo- 
sidereal time. The periods of these modulations 
are: solar diurnal day = 86,400 seconds; sidereal 
day = 86,164.09892 seconds Q; the pseudo- 
sidereal day was assumed to be longer than a 
solar day by the same fraction that a sidereal 
day is shorter than a solar day, or 86,636.54693 
seconds. This pseudo-sdereal wave has zero 
phase at the autumnal equinox, when the side¬ 
real time, the solar time, and the pseudo-sidereal 
time are coincident iQ. In this analysis, the 
pseudo-sidereal time was computed relative to 
the 1988 autumnal equinox, September 22, 1988, 
19/1 29 ™ ux. Second, the live time in each run 
was similarly binned into three histograms. The 
live time for a run was computed as the differ¬ 
ence between the arrival times of the first and 
last muons. The live time was added to the his¬ 
tograms from the time of the first muon to the 
time of the last muon. The rate histograms for 
each run were then computed by dividing the 
contents of the appropriate event histogram by 
the contents of its corresponding live time his¬ 
togram. 

The rate histogram for each run was unpacked 
and the muon rate in each bin, ri, was compared 
to the mean muon rate for that run set, fj, and 
its fluctuation from the mean was computed as 
~ G)/G■ Each Srij for that run was 
then compared with the r.m.s. of the distribu¬ 
tion of fluctuations for the total data set and 
those Srij > (3 x r.m.s.) were cut from the 
analysis. This cut was made to exclude the ef¬ 
fect of the fluctuations found in the long, asym¬ 
metric non-Gaussian tails in the distribution of 
fluctuations for the total data set. These out¬ 
liers, which comprise much less than 1 % of the 
data, mostly result from run starts, run stops, 
sudden data spikes, and other nonphysical ef¬ 
fects. The results of the analysis are relatively 
insensitive 10 %) to this cut over the range 
(2.5 — 5) X r.m.s.. The Srij passing this cut 


three periods. 

At the end of this process, after all runs 
were analyzed, there resulted a set of three his¬ 
tograms with fluctuations from the mean muon 
rate binned in solar diurnal time, sidereal time, 
and pseudo-sidereal time. The contents of each 
bin in these summary histrograms were then 
normalized by the number of entries in that bin. 
The resulting histograms of the normalized fluc¬ 
tuations, AR/Rq, for the solar diurnal, sidereal, 
and pseudo-sidereal periods are shown in Fig¬ 
ures 3, 5, and 6 , respectively. 

IV. PHASE ANALYSIS 

We searched for the first harmonics of the 
solar diurnal, apparent sidereal, and pseudo- 
isidereal modulations by fitting the histograms 
of the fluctuations to the form: 

[Ai?/i?o] = Kcos 2TT[{t - t^ax)/‘24], (13) 

where AR/Rq is the deviation from the average 
muon rate; K is the amplitude of the modu¬ 
lation; t is the solar diurnal, sidereal time, or 
pseudo-sidereal time; and tmax is the time of 
maximum. The results of a x^-fit to the solar 
diurnal and sidereal histograms are given in Ta¬ 
ble 1. The fitted curves are superposed on the 
histograms in Figures 3 and 5. In Table 1, a 
correction of -0.3^ to tmax been applied to the 
sidereal wave fit to account for MAGRO’s ac¬ 
ceptance and overburden, as described in § 11 .A. 
For the pseudo-sidereal histogram, the expected 
pseudo-sidereal modulation from eq.(p^ is given 
in Table 1, with errors computed by adding the 
errors in the yearly and solar diurnal modula¬ 
tions in quadrature. This expected wave has 
been fitted to the pseudo-sidereal histogram and 
the results shown in Figure 6 . 

In addition, Table 1 gives the results of fitting 
the null hypothesis to the three data sets - the 
hypothesis that no modulation is present in the 
deviations from the average muon rate. 

A. Solar diurnal modulation 

Table 1 shows that the solar diurnal modu¬ 
lation is a statistically significant effect. The 
origin of this modulation is most likely the daily 
atmospheric temperature variations at approxi¬ 
mately 20 km, the altitude of the primary cosmic 
ray interactions with the atmosphere. 

We have used the atmospheric temperature 
data from the European Gentre for Medium- 


TT 



12h, and 18h UT, derived using hydrostatic 
equations from satellite observations of the up¬ 
per atmosphere to test the hypothesis that the 
solar diurnal modulation is due to temperature 
effects. We retrieved temperatures in a 25 point 
grid centered in the Gran Sasso region (41.5° N 
- 43.5° N; 12.5° E - 14.5° E) for the 7 atmo¬ 
spheric depths used to calculate Te/f iQ]. The 
result is shown in Fig. 4 superposed on the solar 
diurnal wave. The errors are obtained from the 
dispersion of the hourly data around their aver¬ 
age values for the 7 year period included in the 
analysis. The excellent agreement between the 
observed solar diurnal wave and that expected 
from the daily temperature wave is apparent, 
thus offering strong evidence for the meteoro¬ 
logical origin of the solar diurnal wave in the 
underground muon rate. MACRO is the deep¬ 
est experiment to report this effect. 

B. Sidereal modulation 

Table 1 shows that the sidereal modulation 
is also a statistically significant result. The ob¬ 
served amplitude of K = 8.1 x 10““^ is signif¬ 
icantly larger than expected for the Compton- 
Getting effect due to the solar system motion 
with respect to the LSR, K = 2.6 x 10“"'’; the 
maximum observed phase of t^nax = 23.2^ is also 
significantly different from the expected t^ax = 
17.8^. This suggests that we have found another 
modulation in addition to expected Compton- 
Getting modulation due to solar system motion. 

C. Pseudo-sidereal modulation 

As shown in Table 1, the fit of eq. ( 0 ) to the 
psuedo-sidereal histogram has the same statisti¬ 
cal significance as the null hypothesis. This re¬ 
sult is not unexpected. Since the significance of 
the solar diurnal and sidereal modulations are 
detections at or near the limit of the detector 
statistics, the pseudo-sidereal modulation is too 
small for a statistically significant signal. 

V. RESULTS 

Since the pseudo-sidereal modulation is of 
low statistical significance, the correction to the 
apparent sidereal modulation for the pseudo- 
sidereal modulation is not required. After cor¬ 
recting the sidereal modulation for the solar mo¬ 
tion relative to the LSR, the final value for the 
sidereal modulation of the underground muon 
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AT* = (8.2 ± 2.7) X 10 and a phase of max¬ 
imum of tmax = 0.4 ± 1.1 hr. 

The comparison of the result found in this in¬ 
vestigation with the results of other experiments 
in the primary cosmic ray energy range Ep > 1 
TeV are shown in Fig. 7. In the upper panel 
the amplitudes of the modulations are shown. 
To compare the results of different experiments, 
the reported amplitudes have been corrected for 
the effective declination of the experiment, c.f. 
eq(||) and to the LSR l^j, |^. Included in the 
amplitude is a systematic error of 7% to account 
for the uncertainty in Peff- In the lower panel, 
the sidereal hour of the maximum amplitude of 
the modulation, or the right ascension of the 
maximum signal, is compared with results from 
other experiments. Again, the phase has been 
corrected to the LSR. In this figure, the y-axis 
labels for 18^ to 24^ have been replaced by -6^ 
to 0^. These figures clearly show that the re¬ 
sults found in this investigation are consistent 
with other experimental determinations in the 
energy range 10'^^ — 3 x 10'^"' eV, and that these 
results provide a bridge between underground 
experiments and extensive air shower arrays. 

The interpretation of the sidereal wave is not 
straightforward due to the complicated nature 
of cosmic ray propagation through the galactic 
magnetic field at particle energies in the range 
lO'^^ — lO'^^ eV. However, some understanding 
can come from mapping the observed sidereal 
wave into galactic coordinates. We have realized 
this mapping by first recognizing that MAGRO 
is most sensitive to the declination equal to 
6 ~ Peff ~ 37°, and then tracing out the path 
of this declination in sidereal time. Once this 
path has been determined, we transform it into 
galactic coordinates |^. This path is shown in 
Fig. 8. Onto this path we superpose the sidereal 
wave as follows: positive wave amplitudes are 
shown as crosses and negative wave amplitudes 
are shown as open circles. The position on this 
path of the maximum wave amplitude is shown 
as a filled star. To show the effect of the largest 
possible variation in /3e//, the position of the 
maximum wave amplitude for /3e// = 15°, 64°, 
the outlying peaks in the declination distribu¬ 
tion in Fig. 1, are shown as filled triangles. In 
addition we show the direction toward which the 
sun is moving as a result of differential galactic 
rotation as a filled circle and the direction of the 
Perseus spiral arm as a filled square. This figure 
shows that the wave intensity is positive over a 
wide range in galactic longitude, I = 80° — 150°, 
when MAGRO is pointing towards the galactic 
plane, and is negative when MACRO is point¬ 
ing to high galactic latitudes. This pattern is 



sphere 00- However, it is difficult to ascribe 
the wave to a specific direction and/or “source”, 
as for instance, the direction of galactic rotation 
or the direction of the local magnetic field . 
This is primarily because the wave is at the very 
limit of MACRO statistics. Among the possi¬ 
bilities for the origin of the wave are: particle 
capture in magnetic traps |l^ , a close-by super¬ 
nova [0, ^ ; or the motion of the solar system 
though a spherical distribution of halo cosmic 
rays. 
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FIGURE CAPTIONS 

FIG. 1. Upper panel: Hour angle distribution 
for the muons used in this analysis. The mean 
of the distribution falls at —7.7°. Lower panel: 
The declination distribution for the muons used 
in this analysis. The broad distribution has 
peaks at S = 15°, 42°, and 64°. The mean of 
the distribution falls at 6 = 37°. 

FIG. 2. Percentage deviations of the monthly 
average muon rate from the yearly mean, 
[AR{tq)/R o]yr (solid circles), computed from 
data collected by MACRO during the years 
1991-1997, compared with percentage devia¬ 
tions of the monthly average effective tem¬ 
perature from the yearly mean, AT^ff/ < 
Teff > (open triangles), during the same 
period from balloon data provided by Ispet- 
torato Telecommunicazioni ed Assistenza Volo 
Dell’Aeronautica Italiana. The computation of 
Teff is described in Q. A small offset between 
the monthly means has been introduced for clar¬ 
ity. 

FIG. 3. Deviations of the muon rate from the 
mean muon rate binned according to the local 
solar diurnal time at the Gran Sasso. Super¬ 
posed is the best-fit curve of the form eq. (H 
representing the modulation. 

FIG. 4. Superposed as open triangles onto 
the solar diurnal modulation of the mean muon 
rate, AR/R, shown in Fig. 3 are the deviations 
from the daily mean of the effective tempera¬ 
ture, ATeff/Teff, at 0 , 6 , 12 , and 18 hours. 
The temperature data are taken from the Eu¬ 
ropean Cetre for Medium-Range Temperature 
Eorcasts (ECMWF) on a grid centered on the 
Gran Sasso. 


of the amplitude; Lower panel: measurements of 
the phase maximum of the first harmonic in the 
sidereal variations. Filled symbols represent un¬ 
derground experiments and open symbols repre¬ 
sent EAS experiments. Symbol key: large filled 
circle = MACRO (this paper); inverted filled tri¬ 
angle = Poatina |Q; filled smare = Matsushiro 
||l4[| ; filled triangle = Utah |Q; asterisk = Hong 
Kong Ip^ ; small filled circle = Baksan [ p^ ; filled 
star = Tibet |17t| ; open cross = Baksan EAS 
p^ ; open triangle = Kamioka Q|; open square 
= Norikura p^ ; open diamond = Musala Peak 
]2^ ; open cirele = EAS-TOP Q; open star = 
Linsley 

PIG. 8. The sidereal wave detected by 
MACRO mapped into Galactic coordinates. 
Positive wave amplitudes are shown as crosses; 
negative wave amplitudes are shown as open cir¬ 
cles. The wave’s maximum amplitude, corrected 
for the motion of the sun with respect to the 
Local Standard of Rest is shown as a filled 
star, assuming 6 = 37°. The upper filled tri¬ 
angle shows the maximum wave amplitude as¬ 
suming 6 = 64°; the lower filled triangle shows 
the maximum wave amplitude for S = 15°. The 
direction toward which the sun is moving as a 
result of differential Galactic rotation is shown 
as a filled circle and the direction of the Perseus 
spiral arm is shown as a filled square. 


FIG. 5. Deviations of the muon rate from the 
mean muon rate binned according to the local 
sidereal time at the Gran Sasso. Superposed is 
the best-fit curve of the form eq. &) represent¬ 
ing the modulation. 


FIG. 6 . Deviations of the muon rate from the 
mean muon rate binned according to the local 
pseudo-sidereal time at the Gran Sasso. Super¬ 
posed is the expected pseudo-sidereal modula¬ 
tion from eq.(p^). 


FIG. 7. Summary of sidereal wave searches at 



TABLE I: Results of Histogram Fits 


Period 

K 

tmax (hr) 

xVDoF 

xVDoF 

(null hypothesis) 

Solar Diurnal 

(0.88 ±0.26) X 10“^ 

17.8 ± 1.2 

3.4/5 

14.6/7 

Apparent Sidereal 

(0.82 ±0.27) X 10“® 

23.2+ ± 1.3 

4.6/5 

13.8/7 

Pseudo-sidereal 

(0.0097 ± 0.0029) x 10"® 

4.8 ±2.5 

6.9/7 

6.9/7 


^ corrected for MACRO response in sidereal hour 
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